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We propose a description of dark energy and acceleration of the universe in extended super- 
gravities with de Sitter (dS) solutions. Some of them are related to M-theory with non-compact 
internal spaces. Masses of ultra-light scalars in these models are quantized in units of the Hub- 
ble constant: = n . If dS solution corresponds to a minimum of the effective potential, 
the universe eventually becomes dS space. If dS solution corresponds to a maximum or a saddle 
point, which is the case in all known models based on N = S supergravity, the flat universe 
eventually stops accelerating and collapses to a singularity. We show that in these models, as 
well as in the simplest models of dark energy based on = 1 supergravity, the typical time re- 
maining before the global collapse is comparable to the present age of the universe, t = 0(10""^") 
years. We discuss the possibility of distinguishing between various models and finding our 
destiny using cosmological observations. 
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I. INTRODUCTION 

A. Dark energy and our future 

Recent observations of Type la supernovae and the 
CMB shovif that the Universe is accelerating, and it is 
spatially flat (fitot = + = !)■ Approximately 
0.3 of the total energy density of the universe po ~ 
lO^i^^Mp 10~^^ g/cm'^ consists of ordinary matter 
(JIm ~ 0.3), and 0.7 of the energy density corresponds 
to dark energy {Hd « 0.7), see [|l|,|| and references 
therein. 

One can interpret the dark energy pjj either as the 
vacuum energy (cosmological constant) A ~ 0.7 po, or 
as the slowly changing energy of a scalar field (j) with 
a vacuum-like equation of state pn = w pn, w « — 1. 
In either case our universe is supposed to be acceler- 
ating and rapidly approaching de Sitter (dS) regime. 
Therefore it would be tempting to conclude that our 
universe in the future is going to expand exponen- 
tially for an indefinitely long time, a ^ e^*, even if it 
is closed. 

However, dS regime may be transient, and the fu- 
ture of the universe may be quite different. For exam- 
ple, in most of the models of dark energy it is assumed 
that the cosmological constant is equal to zero, and 
the potential energy V{4>) of the scalar field driving 
the present stage of acceleration, slowly decreases and 
eventually vanishes as the field rolls to (f> = oo, see 
e.g. H]-]^]. In this case, after a transient dS-hke stage, 
the speed of expansion of the universe decreases, and 
the universe reaches Minkowski regime. In both eases 
(dS space and Minkowski space) life may survive for 
an extremely long time until all protons decay and the 
energy resources of the nearby part of the universe are 
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depleted. 

However, there is another possibiUty, which for a 
fong time did not attract much attention. It is quite 
possible that V{(j)) has a minimum at V{(t)) < 0, or 
that it does not have any minimum at all and the 
field (j) is free to fall to V{4>) = — oo. In this case the 
universe eventually collapses, even if it is flat ^-|T^. 

Thus, depending on the choice of the model describ- 
ing dark energy, the Hat universe with = 1 may 
eventually become dS space, or Minkowski space, or 
collapse. (It can never become AdS space with en- 
ergy density dominated by a negative cosmological 



constant 
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B. Dark energy in extended supergravity 

The main goal of this paper is to investigate the 
possibility to describe dS space and dark energy in 
supergravity. We will mainly concentrate on the ex- 
tended gauged supergravities, although we will also 
discuss some models based on A'^ = 1 supergravity. 
The extended supergravities are most interesting for 
two reasons. First of all, they have much closer re- 
lation to M/string theory than the ordinary N=l su- 
pergravity. In particular, the maximal d = 4 N = 8 
supergravity has the same amount of supersymmetries 
as M-theory, 8 x 4 = 32 (there are 32 supersymmetries 
in d = 11 M-theory). Also, in the context of theories 
with extra dimensions the standard N=l supergrav- 
ity (with 4 supercharges) simply does not exist in the 
bulk; the smallest possible supersymmetry in d = 5 
has 8 supercharges, 2x4 = 8 and therefore A^ > 2 in 
d>5. 

In A^ = 1 supergravity it is easy to construct a 
nonvanishing scalar potential (the _F-term potential) 
by choosing a proper superpotential W. The freedom 
of choice is almost unlimited, which does not provide 
us with strict guiding principles that could help us 
to construct a realistic theory. Meanwhile, in the ex- 
tended supergravity one does not have this freedom. 
Nonvanishing potentials appear only after gauging of 
some of the global symmetries of the theory. In a cer- 
tain sense, the total potential is similar to the D-term 
potential in A^ = 1 theory: it disappears for vanishing 
gauge coupling g. 

Recently it was found that one can describe dark 
energy in some d = A extended gauged supergravities 
that have dS solutions ||l^,|l^. Some of these super- 
gravities solve the equations of motion of M-theory 
with non-compact internal spaces [ p^ . 

These dS solutions correspond to the extrema of the 
effective potentials V{4>) for some scalar fields </). An 
interesting and very unusual feature of these scalars 
in all known theories with A" > 2 is that their mass 
squared is quantized in units of the Hubble constant 
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Hq corresponding to dS solutions: ^ = n, where n 



are some integers of the order 1. This property was 
first observed in ||l^ for a large class of extended su- 
pergravities with unstable dS vacua, and confirmed 
and discussed in detail more recently in [p^ with re- 
spect to a new class of A^ = 2 gauged supergravities 
with stable dS vacua pol . The universality of the re- 
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lation %^ = n may be attributed to the fact that ^ 

.° . . ° 

is an eigenvalue of the Casimir operator of dS group. 

Hi- 

The meaning of this result can be explained in the 
following way. Usually the effective potential near its 
extremum can be represented as V{(j)) = A + m'^cjP /2, 
where A and are two free independent parameters. 
However, in extended supergravities with A > one 
always has = nA/3, where n are integers (we are 
using units Mp = 1). Taking into account that in 
ri = A/3, one has, for |0| < 1, ^(0) = 
n02/6) = 3iJ^(l + 7102/6). 



dS space 

A(l -I- n(j)'^/6) = 3H'^{1 + 710^/6). In particular, in 
all known versions of A^ = 8 supergravity dS vacuum 
corresponds to an unstable maximum, 
ll|,ll|], i.e. at 101 < 1 one has 



F(0)=A(l-02)=3i/2(^„^2) 



(1) 



Meanwhile, for the N — 2 gauged supergravity with 
stable dS vacuum found in |20| one has = 6Hq for 
one of the scalars, and for |0| <C 1 one has 



F(0)=A(l + 02)=3i/2(l+02) 



(2) 



In ||l3| , [l^ it was explained that the 'fast-roll' regime 
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with — 0(1) is suitable for the late inflation de- 
scribing the present stage of acceleration of the uni- 
verse. If one takes A ^ po ^ 10^^^° in units Mp = 1 
and Hq ~ 10"^'^ eV, corresponding to the present 
stage of expansion of the universe, then the quantiza- 
tion rule implies that there are ultra-light scalars with 
the mass of the order \m'^\ ^ Hq ^ {10~^^eV)'^. Here 
one should distinguish between the time-dependent 
Hubble constant H^{t) — {pM + Pd)/'^ and its value 
— A/3 in dS regime where pM = and the 
dark energy field (p stays at the minimum/maximun 
of its potential. However, at the present time, with 
« 0.7, one has H{t) = 0{Ho). 
In the early universe the ultra-light scalar fields may 
stay away from the extrema of their potentials; typi- 
cally they 'sit and wait' and start moving only when 
the Hubble constant H{t) determined by cold dark 
matter, decreases and becomes comparable to |m|. We 
will see that this could result in noticeable changes of 
the effective cosmological constant during the last 10 
billion years. The existence of this effect can be veri- 
fied by observational studies of the acceleration of the 
universe. 

Extended supergravity may lead not only to po- 
tentials with dS extrema, but also to exponential po- 
tentials. We will show that a particularly interesting 
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potential V ~ e^*^ can be derived in several models 
based on extended supergravity. Despite the current 
lore 0, the theory with this potential can describe 
dark energy and the present stage of acceleration of 
the universe, and it does not suffer from the problems 
with the existence of the event horizon discussed in 

An important feature of extended gauged super- 
gravities is the fact that quantum corrections to the 
cosmological constant as well as to the ultra-light 
masses are related to the value of the cosmological 
constant A defining the scale of SUSY breaking. For 
A ^ 10~^^'^Mp these quantum corrections are very 
small. Thus extended supergravities provide an ex- 
ample of a model where ultra-light scalars |to^| ~ A 
naturally appear and are protected against quantum 
corrections. 



C. Supergravity and the fate of the universe 

The relation |to^| ~ appears not only in ex- 
tended supergravity. It is rather common in iV = 1 
supergravity for the moduli fields that have vanish- 
ing mass at Mp — > oo [Q. In contrast to the ex- 
tended supergravities, in the = 1 case it is possi- 
ble to avoid this relation: one can take a non-minimal 
Kahler potential, fine-tune the superpotcntial, and/or 
introduce non-trivial D-terms, what will modify the 
mass/Hubble ratio in a significant way. Indeed, in 
the supersymmetry breaking hidden sector one should 
avoid the relation |m^| ^ Hq, in order to avoid huge 
cosmological constant A > (lO'^GeV)'*. However, 
there is no need to make this fine-tuning in the dark 
energy hidden sector^. 

In particular, the relation |m^| ~ occurs in 
the so-called supergravity quintessence model with dS 
minimum . Although the potential of a truncated 
model in this case is simple, the supergravity model 
p4| ] is rather complicated, involving many fields with 
non-minimal Kahler potentials and a set of additional 
assumptions. In order to study a much simpler toy 
model for dark energy in A^ = 1 supergravity we will 
consider a Polonyi-type model of the dark en- 
ergy hidden sector. It has a minimal Kahler potential 
and the simplest superpotcntial W{z) = y?[z + fi). 
The parameter /i should be taken extremely small, 
fj,'^ ^ po ^ lO^^^^Mp. In this case there is no need 

to make the standard fine-tuning (3 = 2— VS in or- 
der to avoid the large cosmological constant. For 
|/3| < 2 — VS the potential has dS minimum with 



* A hidden sector for quintessence models in A'^ = 1 super- 
gravity, different from the supersymmetry breaking hidden 
sector, was proposed in p3,EJ|. 



A - +10^120^,^4 larger values of /3 = 0(1) the 
potential has a minimum with A ~ —\Q~^^^Mp. As 
we will see, in both cases this model, just like the 
models based on A^ = 8 and N = 2 supergravity, can 
describe the present state of acceleration of the uni- 
verse. However, the future of the universe does depend 
on the choice of the model. 

Another interesting model is the axion quintessence 
|2^,|2^. In the M-thcory motivated version of this 
model proposed in one has V{(j)) ~ Acos((/)//); 
for / = 0{Mp) one finds = y'/(o) = ~0{H^). As 
we will show, this version of the axion quintessence 
model can successfully describe the stage of accelera- 
tion of the universe, but, just like the N = 8 models, it 
leads to a global collapse of the universe in the future. 

In order to obtain a fully realistic model of dark 
energy, one would need to construct theories involv- 
ing the observable sector, the hidden sector respon- 
sible for supersymmetry breaking, and the hidden 
dark energy sector. This is a complicated and as 
yet unsolved problem. The supersymmetry breaking 
scale in the dark energy hidden sector is of the or- 
der M^^^y ~ 10^12 GeV . Meanwhile, in the super- 
symmetry breaking hidden sector the corresponding 
scale is M^usy ~ 10^ GeV or even much greater. This 
means that the difference between these two types of 
supersymmetry breaking is more than 15 orders of 
magnitude. Even though the fields from each of these 
sectors may interact with each other only gravitation- 
ally, this interaction may be strong enough to alter 
the important relation m = 0{Hq) for the ultra- light 
scalars. This problem and various ways to address it 
were discussed in p^ , p^ , p^ , p7t . 

In this respect, extended supergravity may be par- 
ticularly interesting as the dark energy hidden sector if 
the mysterious mass quantization rule m? = n Hq has 
some fundamental meaning and remains stable with 
respect to the interaction of the ultra-light scalars 
with the fields from the observable sector. One may 
even argue that the reason for using extended N > 2 
supergravities is due to the nature of gravitational and 
vector fields that may live in five dimensions or higher, 
where A'^ = 2 is the smallest supersymmetry available. 
However, realistic models of supersymmetry breaking 
in the context of supergravity, branes and extra di- 
mensions are yet to be developed. 

For the time being, one may consider the simple 
models of dark energy based on supergravity as toy 
models with some interesting and very unusual fea- 
tures that could be studied by cosmological observa- 
tions. 

In particular, if dS vacuum corresponds to a mini- 
mum of the effective potential, as in the N = 2 model 
of Ref. (20) and in the Polonyi model with < 2-^3, 
the universe asymptotically approaches a stable dS 
regime. 

On the other hand, if the effective potential is neg- 
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ative at the minimum, V{(j)Q) < 0, as in the Polonyi 
model with |/3| > 2 — \/3 and in the axion quintessence 
model with V{(p) ~ Acos(0//), or if it is unbounded 
from below, as in all known versions of = 8 super- 
gravity admitting dS solutions, the flat = 1 universe 
eventually collapses. The typical time remaining be- 
fore the collapse in these models is 0(m~^) ~ Hq^. 
Since the total age of the universe now is also given 
by 0{Hq^), in this class of models the time remaining 
before the global collapse is of the same order as the 
present age of the universe, ^BigCiunch ^ 10^*^ years. 

Thus, in this paper we will study the supergravity 
models which are able to describe nicely the present 
and the past, and may predict an ultimate collapse 
that may occur within the next 10-20 billion years. 
We will show that some of the models predicting grav- 
itational collapse within the next 5 billion years can 
be ruled out by observational data. The investigation 
of the stability of the universe on a greater time scale 
will require a much more detailed observational study 
of the present stage of acceleration of the universe. 

II. SUPERGRAVITIES WITH DE SITTER 
SOLUTIONS AND DARK ENERGY 

In this paper we will consider late stages of evolu- 
tion of universe, as predicted by inflationary theory. 
This means that the observable part of the universe is 
flat and nearly homogeneous, with the metric 

ds^ ^dt^ ~a^{t) dx^ . (3) 

We will assume that in addition to cold dark matter 
with energy density pM the universe also contains a 
slowly changing field 4>. This field may take different 
values in different parts of the universe, but inside our 
horizon the field is almost exactly homogeneous. 

Another assumption is that this field has an ex- 
tremely flat effective potential V{<j)), and the en- 
ergy density of the post-inflationary universe was 
dominated not by this field, but by radiation and, 
later on, by cold dark matter with density pnit) = 

Pa/(^o) (t^) • Therefore even if the field rapidly 
moved in the early universe, eventually it lost its speed 
due to the effective 'friction' related to the redshift of 
its kinetic energy in an expanding universe. Since that 
time the field sits and waits at some point (j)Q until the 
expansion of the universe slows down and the field be- 
comes free to move again. 

The initial position ipo of the field at this stage is 
determined by the processes in the early universe. If 
one assumes that is extremely flat, and that it 

was flat during inflation, then (j)Q can take all possible 
values in different parts of the universe, with prob- 
ability that does not depend on 0O: see e.g. l29|j30|]. 



It is also possible that in the early universe the ef- 
fective potential was quite different. It could have a 
deep minimum at some particular value 4>o, and then 
the position of the minimum could change when the 
density of matter in the universe becomes sufficiently 
small [ p2[ . For simplicity, in this paper we will assume 
that this is not the case and that V{(f>) was always flat. 
However, this is an important and interesting issue to 
be studied in the future. 

Let us consider for a moment an important case 
where V{(j)) = A = const, or assume that ^(0) does 
depend on but 0o corresponds to an extremum of 
the effective potential, V'{(j)Q) = 0. In this case (if one 
ignores quantum fluctuations) the scalar field remains 
at the same point 0o even after the density of cold 
dark matter becomes small. Gradually the density of 
the universe becomes dominated by F(0o) = A. 

At the stage when the universe was cold dark matter 
dominated, its scale factor obeyed the simple equation 
a{t) ^ t^/'^, and the Hubble constant a/ a was given 
by H{t) = ^. Once the universe becomes dominated 
by T^(0o) = A, this regime switches to a{t) — e^"*, 
where Hq = ^A/3 (in units Mp = 1). This is the 
regime of late inflation that can be described in terms 
of dS space. 

De Sitter space is defined as a hypersurface in a 5d 
space -T2 + ^(x^ +y'^ + Z"^ + W'^) = H^^. Here c is 
the speed of light (later we will use the units c = 1). 
In terms of length scale we have 

~ + {X^ + + + W^) ^ Rl . (4) 

Here Rq — cHq is the event horizon of dS universe: 
an observer in an exponentially expanding universe 
a{t) ^ e^°* sees only those events that take place at 
a distance no farther away than Rq. 

Note that this regime begins at pj^j ^ A. Until that 
moment one can use the relation H = After that 
H decreases a bit and becomes constant, H ~ Hq. 
This implies that the switch between the two regimes 
occurs at the cosmological time t ~ 0.5Hq^. Nu- 
merical investigation shows that the present moment, 
when the total density of cold dark matter is about 
0.3/90 and A is about O.Tpo, corresponds to the time 
t « O.SHq^ after the end of inflation. In this regime 
Ho = a/aTS « V0.7po/3 ~ VO^H, where H is the 
present value of the Hubble constant. Therefore one 
can re-express the present age of the universe in terms 
of iJ: i « 0.8H^^ « OMH-^. The last number will 
often appear in our figures since we will measure time 
in units of the present value of H~^. 

The best estimate of the present age of the universe 
is between 13.6 to 14 billion years |]. This yields 
the value Hq^ ~ 11 billion years, which corresponds 
to the size of the event horizon in dS space Rq = 
cH^^ w 10^8 cm. This corresponds to IQ^^lp which 
is the minimal size of dS hyperboloid, the size of its 
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throat. Here L 



10 



-33^ 



"cm. The corresponding scale 
of energies is IQ-^^Mp - lO'^^ eV. 

Now let us see whether this cosmological picture can 
be related to supergravity. A generic four-dimensional 
gauged supergravity which will be used for dark en- 
ergy hidden sector has part which includes gravity 
coupled to 2n scalar fields and a potential 
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Here Gy {4>) is the metric in the moduli space and we 
have not included vector fields and fermions. 

As explained in ||T^, in all known gauged extended 
supergravities possessing dS solutions, the eigenval- 
ues V" of the mass matrix (m^)^ in dS extremum are 
proportional to the value of the potential V in the 
extremum. This corresponds to integer values of the 
ratio of eigenvalues of the mass matrix to Hq in the 
range between 12 and —6. The mass matrix at the dS 
critical point is defined as follows: 



~iik 



(5) 



Here (/)*^ are the critical values of the scalar fields cor- 
responding to the extrema of V . 

We assume that there is a dark energy Lagrangian 
given by extended supergravity with dS solution and 
that in addition to the scalars <jf representing the dark 
energy of the universe there is also the usual cold dark 
matter energy density p^v/ contribution on the right 
hand side of the Friedmann equations: 



3^. 
a 



dV 







Hit) ^ - = 
a 



PM + V + 



kin 



(6) 



(7) 



Here a{t) is the scale factor of the flat FRW met- 
The kinetic part of the dark energy, 



ric, PM 



Ekin{4>), is given by Ekm = lGtj{<p)(f)'(j>> and r*^((/)) 
is the ChristofFel symbol in the moduli space defined 
by the metric Gij {(f)) . 

In the linear approximation for the fields <j) which 
are close to their critical values (t)cr (i-e. to the ex- 
trema of V{(l))) the scalar field equations are reduced 
to -I- m5'4,^ + {m'^)'jS(l>> = where dcj)' = - 4>i^. 
This shows the universal character of dS critical points 
in all known gauged supergravity models: all scalar 
equations near dS critical point take the form 



-I- 3HoS(f,' + nH^S(f>' = 



(8) 



where n — jri, the Casimir operator of dS group, 
takes values from —6 up to -1-12. 

The energy density po and the pressure Pd for the 
supergravity dark energy are given by 



PD = Ek^n +V= -Gy (<^)(/.>^' + y (</)) , (9) 



Pd - Ekrn -V= -G„- - 1/(0) . (10) 



The total energy includes also the energy of matter 



PT = PM + Ekin + V . 



(11) 



A dark energy (matter) density is given by a ratio of 
the dark energy (matter) to the total energy. 



PD_ 

PT 



M 



PM 
PT 



(12) 



Experimentally now ^Im ^ 0.3 and iljj ^ 0.7 and 
= SIm + f^D '^-^ 1. Thus when we plot the time 
dependence of ^Id the present time will be specified 
by — 0.7. Another important characteristics of 
the dark energy is its pressure-to-energy ratio defining 
dark energy equation of state Pd = wdPd'- 



WD = 



Ekin — V 



PD 



E, 



kill 



V 



(13) 



In what follows we will skip the index Dinwu and use 
w for dark energy wd, since we consider only recent 
period when radiation is irrelevant and pm — wm = 0. 

The parameter w in general strongly depends on 
time t. It is convenient to introduce an average value 
of w which is useful for understanding an integrated 
influence of w on observations: 



J da flDia) w{a) 
J da flnia) 



(14) 



Note that the deceleration parameter (by the time 
when radiation is negligible) is given by 



-[nM + noii + ^w) 

\ ilnW . 

2 2 



(15) 



Here we used the fact that wm = 0. At present Qq ~ 
i -I- w is close to —0.5 and therefore w is close to —1. 



A. Collapsing universe 

In all models with dS maximum or saddle point 
the universe will eventually collapse, i.e. the scale 
factor will shrink to zero. At some point in time t the 
potential becomes negative and cancels the positive 
contribution of pM and Ekin in the expression for the 
Hubble parameter: 



H\i) = 



PM + V + Eki 



= . 



(16) 



At this moment the universe stops expanding; the 
Hubble parameter becomes negative for t > t, which 
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inevitably leads to a collapse |p|-|l7[. In case of col- 
lapsing universe, we have to use an alternative form 
of one of the Friedmann equations. Instead of Eq. 
we will use 



a _ -pM + V - 2Ekr. 
a 3 



(17) 



since Eq. (|7|) is not well defined in the case where 
H vanishes and becomes negative. The equation of 
state function w(t) will grow from the initial value —1 
to infinity somewhere in the future, well before the 
collapse, and then, near the singularity, it usually ap- 
proaches the value -1-1 corresponding to the stiff equa- 
tion of state p = p |16 1^. 



^initial ^ 



today _ 



= -1 



future 



-1. 



B. Future dS universe 

In all models based on extended supergravities with 
dS minimum of the potential, w will start at —1 at the 
time when all ultra light fields are frozen somewhere 
at the slope of the potential, so that = ^- 

the time when the Hubble parameter decreases to the 
value comparable with the mass of these fields, they 
will start moving towards the minimum of the poten- 
tial with some velocity. This will raise the value of 
w above —1. However, eventually all fields will ap- 
proach the minimum and stabilize at E'^ln"''^^"^ = 0, 
perhaps after some oscillations. Again, the attractor 
value w — ^1 will be reached in the long term fu- 
ture. Depending on initial conditions of the fields, at 
present w may be positive, or negative or vanishing 
and w takes any value above —1. This would mean 
that during the last Hubble time it could be increasing 
monotonically or first increasing and later decreasing. 



^initial ^ 



today _ 



-l + 5a 



future 



-1. 



For general supergravity models, we find numerical 
solutions of the Friedmann equations for the universe 
evolution during the last Hubble time. We pay partic- 
ular attention to time dependence of ^.^{t) and w{t) 
since they can help to confirm or rule out the class of 
dark energy models with future dS space or collapse. 
All models are adjusted so that today corresponds to 
Iq = 0, t is negative in the past and positive in the 
future. We use units in which = 1 at present and 
assume that = 0.7, = 0.3 today. 

For each model we also present the evolution as a 



function of the redshift z where 



a(to) 
a(t) 



\ + z and 



Ztoday = 0. We Only consider non-negative z > 0, 
which correspond to the past and present since only 
such z are subject to cosmological observation. 



III. GAUGED SUPERGRAVITIES AND 
COLLAPSING UNIVERSE 

A. TV = 8 supergravity and M-theory 

Our particular interest in the cosmological aspects 
of the theory with maximal supersymmetry is due to 
the common trend to have maximal amount of su- 
persymmetry allowed by experimental and theoretical 
considerations. 

M-theory supergravity or supergravities related to 
low energy string theory, compactified on spaces other 
than torus, may lead to non-trivial potentials. How- 
ever, if the compactification volume is finite the rel- 
evant four dimensional supergravities do not have dS 
solutions 

Still some four-dimcnsional gauged N = 8 super- 
gravities with 32 supersymmetries have dS solutions, 
and some of them are directly related to M/string the- 
ory |l^ . In agreement with the no-go theorems of 
they correspond to solutions of M/string-theory with 
non-compact internal seven (six) dimensional space.|^ 

Recently the meaning of such M/string theories 
with non-compact internal space was reconsidered 
HJ] . The relation between states of the higher dimen- 
sional and the four dimensional theory in such back- 
grounds is complicated since the standard Kaluza- 
Klein procedure is no longer valid. It is nevertheless 
true that the class of four-dimensional dS supergravi- 
ties which we will consider below as dark energy candi- 
dates has a direct link to M/string theory. They solve 
equations of motion of the eleven-dimensional super- 
gravity with 32 supersymmetries. Moreover, these 
theories are perfectly consistent from the point of view 
of the four-dimensional theory: all kinetic terms for 
scalars and vectors are positive definite. 

One of the simplest solutions of lid supergrav- 
ity is given by a warped product space of the four- 
dimensional dS space and the seven-dimensional hy- 
perboloid H'P''^. Our fiducial model with all scalars 
constant is defined by the dS surface (^, and by the 
internal space H^''^ given by the equation tjabz^z^ = 
a i?Q in an eight-dimensional space. Here the constant 
a depends on p, q. The metric r]AB is constant and 
has p positive eigenvalues and q negative eigenvalues, 
p + q = 8. 



^ There exists only one known version of quintessence 
derived from M-theory Q. It starts with 50(3,3) 
gauged supergravity in five dimensions which has dS solu- 
tion, related to eleven-dimensional supergravity with non- 
compact six-dimensional space. The five-dimensional dS 
solution in this model is unstable so it is likely 

to lead to a universe collapsing in the future, just as all 
N = 8 models considered in this section. 
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The simplest (and typical) representative of dS su- 
pergravities originated in M-theory, has the potential 
[ |3^ with dS maximum: 

y = A(2-cosh\/2</.) (18) 

(for canonically normalized fields 0). 

In this model the extremum of the potential is at 
= 0. V{0) = A, y"(0) = -2^(0), which corre- 
sponds to = —QHq. This is the four-dimensional 
N = 8 supergravity with p = q ~ 4, which has gauged 
SO{A, 4) non-compact group. At the dS vacuum it is 
broken down to its compact subgroup, 50(4) x 50(4). 
The value of the cosmological constant A here is re- 
lated to the value of the gauge coupling g as follows: 

A = 25^ 

Near the extremum the potential can be represented 
by the simple quadratic expression V{<j)) — A(l — clP), 
see Eq. (^. It is unbounded from below, so the theory 
is unstable, and the first idea would be to discard such 
models altogether. However, the potential remains 
positive for |0| < 1, and for small A the time for the 
instability to develop can be much greater than the 
present age of the universe, which is quite sufficient 
for our purposes jl^. In fact, we will see that this 
instability allows us to avoid the standard fine-tuning 
problem plaguing most of the versions of the theory 
of quintessence Q . 

The gauge coupling as well as the cosmological con- 
stant in four-dimensional supergravity have the same 
origin in M-theory: they come from the flux of an 
antisymmetric tensor gauge field strength The 
corresponding 4-form -F0123 in eleven-dimensional su- 
pergravity is proportional to the volume form of the 
dSi space: 

F0123 - \/A Vbl23 • (19) 

According to this model, the small value of the cos- 
mological constant is due to the 4-form flux which 
has the inverse time-scale of the order of the age of 
the universe, H^^ ^ to ^ 10^" years. 

Somewhat more complicated potentials are found 
in the case with p — 5,q — 3 with 5*0(5, 3) gauging 
broken down to 50(5) x 50(3) at dS vacuum, and 
in the case with p = 3,q = 5 with 50(3, 5) gauging 
broken down to 50(3) x 50(5). 

Near dS point all known potentials of the gauged 
N = 8 supergravity have the following universal fea- 
tures ll^ : the ratio ^ of the eigenvalues of the mass 
matrix to the potential at the extremum takes values 
-2, 4/3, 4. This imphes that ^ = {-6,4,12}. In afl 

models there is a tachyon direction with = ~QHq, 
and for the second scalar one finds either = AHq 
or = 12^^0' i-'^- the potential has a saddle point. 



So far no gauged N = 8 supergravities with 32 super- 
symmetries were found with dS minimum. Therefore 
the model with the potential ( p^ ) represents also all 
other models, since near the critical point the tachyon 
has the same mass. 



B. Cosmology in = 8 models 

Since all dS critical points linked to M-theory are 
unstable and the potentials are unbounded from be- 
low, one could expect that these models cannot de- 
scribe the past and current evolution of the universe 
and play the role of the dark energy hidden sector. 
However, in ||l^, |l4j it was found that if A ^ po ~ 
the universe may stay now and during the last 
few billions of years near the top of the scalar potential 
while the scalar field was slowly drifting away from the 
critical point. Here we will study this model in detail 
and solve the Friedmann equations numerically. 

To avoid misunderstandings, we should emphasize 
that the relations = uHq have been derived for 
dS extrema with the energy density fully supported by 
the scalar field potential ¥{(/)): m? = nA/3. In our 
calculations we will assume for simplicity that this 
relation remains preserved in the presence of matter 
fields with the energy density pM- 

In our calculations we will also assume that initially 
the cold dark matter energy density was many orders 
of magnitude greater than the energy density of the 
scalar field 0, PM{to) ^ Pd- In such cases the field 
freezes very rapidly, so we can simply start our cal- 
culations assuming that its initial velocity vanished. 
We considered 4 different initial values of the field </>: 
00 = 0, 0.2, 0.3, 0.35 in the units Mpi = 1. The case 
with 00 = Oj where the field starts and remains at 
the top, corresponds to our fiducial model. The cor- 
responding curves are shown as thick red. The blue 
dashed curves correspond to 0o — 0.35. 

For each of the values 0o = 0, 0.2, 0.3 we numeri- 
cally found different values of A such that for each of 
these three cases one has = 0.7 and iJ = 1 to- 
day. Here H is normalized in units of the presently 
observed Hubble constant. However, we could not do 
it for 00 — 0.35 because the largest value of il^) in this 
case is about 0.65; we presented this case for H = 1 
and — 0.65 for comparison. We plot in Figure 1 
the evolution of the scale factor a(i), in Fig. 2 the 
evolution of Vlo{z), and in Fig. 3 the evolution of 
the equation of state factor w{z). The point t — Q 
{z — 0) corresponds to the present. In all cases ex- 
cept 00 — 0.35 the Big Bang occurred at the moment 
to — 0.95i/^^ in terms of the present value of the 
Hubble constant H ~ 1. If one assumes, according to 
1^], that 1^0 1 ~ 13.6 billion years, one finds that each 
interval At = 1 in our figures approximately corre- 
sponds to 14 billion years. 
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FIG. 1. Scale factor a{t) in the model based on A'^ = 8 
supergravity. The upper (red) curve corresponds to the 
model with (po = 0. In this case the universe can stay at 
the top of the effective potential for an extremely long 
time, until it becomes destabilized by quantum effects 
fl^ . The curves below it correspond to <j)o = 0.2 and 
cpo = 0.3. The blue dashed curve corresponds to 4>o = 0.35. 
The present moment is t = 0. Time is given in units of 
H~^{t = 0) ^ 14 billion years. 

As expected, all 3 models with = 0.2, 0.3, 0.35 
lead to a collapse. Still all results for the cases 
00 = 0.2, 0.3 are compatible with observations. The 
case (/)o — 0.35 is close to be ruled out on the basis 
of current data with Qo = 0.7: in this model it is 
impossible to reach Ojj = 0.7, the largest value of ilo 
in this case is about 0.65. Among all of our models, 
this was the case with nearest collapse in the future 
(7 billion years from now). We cannot claim that this 
case is definitely ruled out due to the cosmological 
observations suggesting that = 0.7 today. How- 
ever, the present value of w for this model is very 
high, which is an additional evidence suggesting that 
this regime is disfavored by observations. For greater 
values of 0o the maximal value of^u becomes consid- 
erably smaller. For example, for (f>Q > 0.4 the maximal 
value of flD is 0.56. This means that the models with 
(po > 0.4, predicting even earlier collapse of the uni- 
verse, contradict observational data. 

On the other hand, the models with (po — 0, (j>o — 
0.2, and (po = 0.3 are consistent with the present ob- 
servational data. We conclude that the N — 8 models 
with the potential unbounded from below are viable 
candidates for the description of the present stage of 
acceleration of the universe if the scalar field 4> initially 
was in the range < ^ 0.3 in Planck units. 

Other models of TV = 8 gauged four-dimensional 
supergravity linked to M-theory ||l9] have the same 
important property that = —QH^ fl^]- Therefore 
for the initial values of fields not far from the top of 
the potential, an analogous behavior is expected. 
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FIG. 2. Dark energy Qd ss a function of redshift z for 
(po ~ 0, 0.2, 0.3, 0.35. The lower (red) curve corresponds 
to the model with (j>o = 0. The present time corresponds 
to 2 = 0. Note that for (po = 0.35 the maximal value of 
0,0 is about 0.65. 
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FIG. 3. Equation of state function w as a, function of 
redshift z for (po = 0, 0.2, 0.3, 0.35. For (po = 0.35 this 
function sharply rises near z — 0. For (po = 0.3 the maxi- 
mal value of w is about —0.65. This could also seem too 
high, but the average value of w in the important interval 
2 2 is below —0.9. For (po — 0.2 the maximal value of 
w is below —0.9. The red line w — ~1 corresponds to the 
model with (po = 0. 

One could argue that fixing A ^ 10^^^" and (j> < 0.3 
in Planck units represents fine-tuning. Of course, fine- 
tuning of such type is present in most of the models 
of dark energy. But can we do better than that and 
avoid fine-tuning altogether? 

Note that one cannot simply add a nonvanishing 
cosmological constant to N=8 supergravity keeping 
the curvature of the effective potential intact. This 
would violate the structure of the effective potential 
V = A(2 - coshV20) and the relation = -2A. 
Once we increase A, the tachyonic mass increases as 
well, and the total lifetime of the universe becomes 
unacceptably small. Indeed, as we have found, the 
total lifetime of the universe in such models typically 
is 0{H-^) - 0(A-i/2). It can be shown that this 
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lifetime is proportional to log (/)q , so it can be made 
much greater than H~^, but only if the initial value 
(po of the field (j) was exponentially close to zero p^ . 
When the field (po becomes very close to 0o = 0, the 
lifetime of the universe becomes determined by quan- 
tum fluctuations. The total lifetime of the universe 
cannot exceed t ^ lO^Hf^^ ^ 10^ billion years even 
for very small (/)o (is) . 

Let us assume that for the development of life of our 
type the universe must have lifetime ^tot ft 14 billion 
years, which is the present age of the universe. This 
would imply that we cannot live in a universe with 
A ^ lO"""^^" since such a universe will be short-living, 
with the total age much smaller than 14 billion years, 
unless one assumes that in the very beginning (/jq was 
exponentially small. 

According to Eq. (|l^), the value of VA is fixed 
by the 4-form flux, ^0123 ~ ■s/AVbi23. In general one 
may imagine that this flux (and, consequently. A), 
as well as the field (f>o, may take different values ei- 
ther in different regions of the universe, or in different 
quantum states of the universe. Consider all combi- 
nations of (f>Q and A compatible with the total lifetime 
of the universe itot > 14 billion years. If one makes 
an assumption that all such combinations are possi- 
ble and equally probable (or equally natural), then 
our investigation indicates that the value of A should 
very close to po ~ 10~^^°, and the most probable time 
before the future collapse of the universe is 0(10) bil- 
lion years |^6| . Depending on various assumptions, 
one finds also that the probability to have flo in the 
interval between 0.5 and 0.9 at the time t ^ lA bil- 
lion years after the big bang is about 15% [^. This 
number may change if we make different assumptions 
about the probability measure, but many reasonable 
assumptions leads to the same conclusion: Finding 
rio in the interval between 0.5 and 0.9 does not re- 
quire the fine-tuning that is the usual problem of many 
models of quintessence. This may provide a solution 
to the fine-tuning/coincidence problem for A and Qjj 
and, simultaneously, a prediction for the typical time- 
scale of the global collapse of the universe in the mod- 
els based on = 8 supergravity. 

C. Collapsing universe in A*' = 4, A'^ = 2 gauged 
supergravities 

The very first gauged supergravity with dS solu- 
tions, breaking supersymmetry spontaneously, was 
found in |3^. The relevant part of the action of 
SU{2) X SU{2) gauged supergravity is 



V" /V = -2 , ^ = -6. This means that the cosmo- 
logical evolution in this model for initial values of the 
fields not far from the dS critical point at \W\ =0 will 
be the same as in the previous case, explained above 
for TV = 8 theory. 

As an example of unstable dS solutions consider 
here N = 2 SU{2) gauged supergravity [|8|. The rel- 
evant part of the action depends on the scalar fields 
i?^ and R^, A = 1,2,3: 

g-^'^L =\r- SAEdf^R^d^R^ + 2g^[3 - 2R^R^Sab 

~ 5AB{eACDR'^R''){eBEFR^R^)] ■ 

This action has two dS critical points. One is SU (2) 
symmetric and Rf^ = 0, Vcr = 37?o = 65^. The 
tachyonic mass is equal to ~Ag^ and = , jp- = 

—3. At the second critical point {R^R^ 5AB)cr = 1, 
{R^R^5AB)cr = 0, Vcr = 3iJ^ = 4g2. There is a 
tachyonic mass squared equal to —85^ and a positive 

mass squared 1617^ and some Goldstone modes. Thus 

2 

we find in the second minimum that ^ = —6, 12. 
Thus we have shown in examples above that despite 

a large variety of potentials in extended gauged super- 

2 

gravity, we always confirm the pattern that ^ is of 

the order one and there are tachyons. The develop- 
ment of instability and all basic features which were 
found by numerical solutions of the Friedmann equa- 
tions in A^ = 8 case will take place in other models 
with unstable dS solution. 



IV. GAUGED SUPERGRAVITIES WITH DE 
SITTER ATTRACTOR 

The models with dS minimum always define ulti- 
mate future dS universe since scalar fields eventually 
reach their constant fixed points at which the poten- 
tial has a minimum, and space-time is a dS-type. The 
simplest (and typical) representative of a potential in 
such theories (for canonically normalized fields <j)) is 

1/ = A cosh V2(t> . (20) 

Here A is the cosmological constant, the value of the 
potential at the minimum. It is related to the attrac- 
tor value of the Hubble constant, A = 3i?Q. At the 
minimum where y = 0, = and V{0) = A > 0. 

A. Axion-dilaton dark energy 



^ 2 



1 



rd,Wd^W 



1- liyP " " 1- |iy|2 

Here the potential depends only on |T4^| and again 



N — 2 supergravities with stable dS vacua con- 
structed in poll have special features: all moduli are 
coordinates of the special Kahler and quaternionic 
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manifold 5T[2, n] iJQ[m]. The ungauged super- 
gravity of this type appear, e. g. in the large ra- 
dius limit of superstring compactifications | ]39| . Upon 
consistent gauging of the non-compact groups directly 
in d = 4, these theories acquire a potential and be- 
come rather complicated, in general. The simplest 
version with 3 vector multiplets has the axion-dilaton 



S'-moduli which parametrize an 



SU{1,1) 
U(l) 



part of the 



manifold and the so-called T and U moduli (Calabi- 
Visentini coordinates, yo,yi ), which parametrize the 
^'-^^^'^^ — part of the special Kahler manifold. 

to focus on the simplest 



SO(2)xSO{2) 

It was suggested in 
part of all these models which consists of axion-dilaton 
only. One reason for that is the simplicity, the other 
reason is that axion-dilaton pair in 4d theory, com- 
ing from M/string theory, have well known global 
SL{2, Z)-symmetry. It includes the Peccei-Quinn shift 
of the axion. The scale of violation of this symmetry 
is related to the scale of the cosmological constant and 
mass of these fields. It is therefore interesting to study 
the evolution of the universe during the last Hubble 
time and the future evolution towards dS attractor 
and see what kind of features of dark energy and ac- 
celerated universe may be described by this model. 
We will also study a more general model with dS fu- 
ture universe which include Calabi-Visentini scalars in 
addition to axion-dilaton fields. 

Axion-dilaton action (without a potential, before 
gauging) in four dimension obtained by dimensional 
reduction from d = 10 is given by 



-1/2 r _ 9'"'d^Sd,S 
(2ImS)2 ■ 



(21) 



Here the real part of the complex modular parame- 
ter S = \f2A - ie-^'t' is the axion field A and the 
imaginary part has the dilaton (j). 

The action for canonically normalized kinetic terms 
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-1/2. _ 



(22) 



is easily obtained either by compactification on a six- 
torus of the heterotic string theory in ten dimensions 
or by a compactification on a seven-torus from eleven 
dimensional M-theory supergravity action. For exam- 
ple, in string theory (j) is related to the dilaton of the 
string theory and A is related to the second rank anti- 
symmetric tensor iJ-field of string theory, after a du- 
ality transformation. These are so-called fundamen- 
tal dilaton and axion. In M-theory is related to 
component of the metric and A is related to the 
third rank antisymmetric tensor field C. The details 
of(i=10,d = 4 connection are described in |^0|. The 
axion-dilaton system in string theory and ungauged 
supergravity has a non-compact SL{2, R) or SL{2, Z) 
symmetry: the actions (|2l|), (|22|) are invariant under 



linear fractional transformations of the modular pa- 
rameter S ^It^. This includes the shift of the 

cb+d. 

axion field by a constant, A ^ A + const. 

In case of interest a consistent gauging directly in 
cZ = 4 produces a following potential of the axion- 
dilaton action lEoll, H: 



V = A [cosh 72(0 - (/jcr) + e^('^+'^-)(A - AcrY' 



(23) 



Here A is the cosmological constant, the value of 
the potential at the minimum. It is related to the 
attractor value of the Hubble constant, A 
The attractor values of the axion-dilaton field S, 



3Hl 



V2Ac 



'V2<f>ar 



At the minimum where V — 0, 



A^Ar. 



V{A„Acr)^A>{). (24) 



Near the dS vacuum solution with the axion-dilaton 
field fixed at the attractor values S = Scr, the masses 
of the dilaton and axion fields are both equal to = 
TO^ = 6 f = 2A. 

Note that the global axion shift symmetry A 
A -\- const is broken by the potential with the fixed 
value of Acr- 

The original parameters used in [ pO| are some par- 
ticular combinations of gauge coupling eo, Fayet- 
Iliopoulos term ei and magnetic rotation angles 9, 
such that 



A = eoci sin 9 , Acr — cot 9/ V2 e 



eo 



ei smf 



(25) 



Here the non-compact 5*0(2, 1) x U{1) group is 
gauged: eo is a gauge coupling for 50(2, 1) and ei 
is a FI term for U{1). The theory has 4 vector 
fields: a graviphoton and 3 vector fields for each of the 
STU moduli. At the dS minimum the non-compact 
50(2, 1) gauge group is broken spontaneously down 
to 5*0(2): the gauge vectors associated with the non- 
compact generators of the 50(2, 1) become massive 
by eating one of the complex moduli. 



B. S,T,U moduli dark energy 

An un-truncated version of the first model in pO| 
has S,yo,yi fields or, equivalently S,T,U fields. We 
will choose a gauge where the Goldstone scalar field 
j/o is vanishing. This will give us additional kinetic 
terms in the action when yi — ere*'' are present: 



(l-cr2)2 

The total potential is 



[{daf+a^idSY 



(26) 
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y = A 



1 + ^^ 

[1 - cr2]2 



cosh V2(0 - 



(27) 



If we put the axion-dilaton fields in the minimum of 
the potential, i.e. at = (j)cr and A — A^r, the re- 
maining action for the a field is 



1 + 



This can be rewritten as 



(28) 



(29) 



This model is a typical representative of iV = 2 super- 
gravity with future dS space. 

We study numerical solutions of the Friedmann 
equations with the action given in (|3^). Our choice 
of initial conditions is (^o = 0, 0.6, 1.0, 10 and in all 
cases 00 = 0. As in previous models of collapsing uni- 
verse, we choose for 'today' = 1 and ^1^1= 0.7. 
We plot in Figure ^ the evolution of the scale factor 
a{t), in Fig. ^ the evolution of Qd{z) and in Fig. 
H the evolution of the equation of state factor w{z). 
The point z = as in collapsing models corresponds 
to 'today'. The red curves (the upper curve in Fig. ^ 
and the lower lines in Figs. ^ and y) correspond to a 
fiducial model where the scalar field is already at the 
minimum of the potential. 



which gives a very simple gauged supergravity moti- 
vated model of the field with the mass-Hubble ratio 
m2 = A = 3Hl 

The model with 3 scalars, 0, A, a is 



1-1/2 r _ 



1 



1 



L = -^R+^[{dcl,r+e'^ndA) 



(l-a2)2 

2^2 
[1 - a2]2 



1 + 



■cosh%/2(0- 0cr) (30) 



sinh V2{<f) - (j)cr) + e^('^+'^'=-) {A - A^rf 




We solved the Friedmann equations numerically for 
the general model with all 3 moduli fields, correspond- 
ing to the action given in eq. (|30|). The solutions of 
these equations for various initial conditions for the 
scalar fields show the attractor behavior of the total 
solution: the scalars eventually reach their attractor 
values, defined by the minimum of the potential, and 
the universe asymptotically becomes dS space. Rela- 
tive importance of various fields during this process 
depends on their initial values. However, we have 
found that if we begin with the values of the fields 
far from the minimum of the effective potential, then 
for a rather broad range of initial conditions the dila- 
ton contribution becomes dominant at the late stages 
of the process. Therefore we will concentrate now on 
the simplest case of the dilaton-dominated universe. 



C. Generic future de Sitter universe 

Let us assume that the (T-field and the axion have 
already reached their fixed points, cr = and A = Acr 
and only the dilaton is not yet there. We also choose 
the vanishing value of 4>cr = for simplicity: 



-1 -0.5 0.5 1 I 

FIG. 4. Scale factor a{t) in the model based on 
N = 2 supergravity with a stable dS minimum for 
(po = 0, 0.6, 1.0, 10. The upper (red) line corresponds 
to <j!>o — 0, the lower (blue) line corresponds to (j>o — 10. 
Blue lines in this figure and in the next two figures practi- 
cally coincide with the corresponding lines for the purely 



1-1/2 r _ 



L = - 



-R+ -{d(t)f - A [cosh%/2(/)] 



(31) 



exponential potential V ^ e 



In each case we find the universe that ultimately 
becomes dS space, forever expanding, never collaps- 
ing. In the past for f < as well as in the future for 
t > for all cases we find almost the same curves for 
a{t) and floit). All models are very close to the fidu- 
cial model with regard to a{t) and riD(i). Only the 
evolution of the equation of state factor w{t) varies 
from case to case. The importance of this fact is that 
it may at least in principle become observable. For 
all models with 0o > 0.6 we find w > at t = 0, for 
(j)Q — 0.6 ?i; = at t = and for (po < 0.6 zi; < at 
t = 0. 

We should note that this model requires fine tuning 
of initial conditions. We began the calculations at 
some (j) = (j)o, and adjusted the value of the effective 
potential in its minimum in order to obtain = 0.7 
at the time ~ 14 billion years. However, this is 
the standard fine-tuning that is present in most of 
the models of quintessence; see a discussion below. In 
this respect it is interesting that the models with dS 
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maximum considered in the previous section may not 
require fine-tuning |3q|. 
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FIG. 5. Dark energy Q,d{z) for (jyo = 0, 0.6, 1.0, 10. 
The lower (red) line corresponds to (j>o = 0, the upper 
(blue) hne corresponds to <^o = 10. 
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FIG. 6. Equation of state w{z) for = 0, 0.6, 1.0, 10. 

The lower (red) line corresponds to 0o = 0, the upper 
(blue) line corresponds to cf>o — 10. 



V. M-THEORY AND DARK ENERGY WITH 
EXPONENTIAL POTENTIALS 

The results obtained above are closely related to in- 
vestigation of dark energy in the theories with purely 
exponential potentials. Indeed, the blue line corre- 
sponding to the initial condition (/jq — 10 (the lower 
line in Fig. ^ and the upper line in Figs. |^ and ^ 
practically does not change when 0o increases even 
further. This result has an interesting interpretation. 
In the limit 0o ^ 1 (for a = and A = Acr) the 
dilaton potential ( p3| ) can be represented as a simple 
exponent: 



V = Ae^^ 



(32) 



In this case increase of 4>q to 0q can be absorbed by 
rescaling A — > A exp V2('/'o ~ 0o)- This is what we 
are doing in our calculations when we are normalizing 



the present Hubble constant H to H ~ 1. As a re- 
sult, for all (f>o in the exponential potential one has the 
same curves n{z) and w{z). These curves practically 
coincide with the blue lines in Figs. ^ and 

This conclusion has important implications. The 
first models of dark energy were based on investiga- 
tion of the theories with exponential potentials Q. 
For V = Ke^'^ different possibilities were studied, de- 
pending on the value of A. In the presence of cold dark 
matter with equation of state p = 0, the interesting 
regions are: 

• For < 3 the system has a late-time attractor 
with fixed point values fi^**'" = 1 and w''**'' = 
AV3-1. 

• For A^ > 3 there is a stable solution of a scaling 
type, a global attractor with, fixed values fi^**'' = 



3/A^ and 



= 0. 



In the first case one does not have il^i = 0.7 in the 
attractor regime. In the second case the power-low 
regime with VLjj Ki 0.7 is possible. However, since VLjj 
is approximately constant, it is difficult to reconcile 
this regime with the big bang nucleosynthesis which 
requires that at the stage of the nucleosynthesis 17 
must be very small |^ . Moreover, the limit w — > is 
not good for quintessence. Therefore the general lore 
is that the theories with exponential potentials cannot 
describe dark energy. 
In our paper (see also 
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42 1) we study a regime 



that occurs before the late time attractor is reached, 
for A^ < 3. We simply check whether we can find such 
parameters of the theory that would allow us to con- 
sistently describe the observational data. And the an- 
swer to this question is positive. If one chooses proper 
initial conditions (the universe dominated by CMB), 
initial values of become negligible, so one has no 
problems with the big bang nucleosynthesis. Mean- 
while at some later stage Vlu can grow and reach the 
desirable value 0.7. Therefore exponential potentials 
can describe the present stage of acceleration of the 
universe. This description requires some fine-tuning, 
but it is not worse than it is in many other models 
of dark energy Our conclusion agrees with the one 
obtained in p2| . 

Another problem associated with exponential po- 
tentials is that they must be sufficiently fiat. It is well 
known that the universe, in the absence of ordinary 
matter, will accelerate only for A^ < 2; for A^ = 2 one 
would have a linear regime a ^ t with a = 0. However, 
in order to explain observational data we do not have 
to assume that the universe will keep accelerating in 
the distant future, so the requirement A < is not 
necessary. On the other hand, our numerical studies 
show that for A > 2.2 the value oi ^Ijj never reaches 
0.7. For A > 1.7 one can have flo = 0.7, but one has 
w > —0.6 at present, which is perhaps too high B3]. 
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We may conclude therefore that for < 3, when 
the system has a late-time attractor, the reason- 
able description of the present and past can be eas- 
ily achieved at the stage before the late attractor is 
reached. Eventually, in the long term future the 
fixed point values will be reached, flo 1 and 
m;^AV3-1. 

Recently it was pointed out that it might be dif- 
ficult to describe eternally accelerating universe with 
< 2 in terms of M/string theory because it has 
event horizon |2^. In case A^ = 2 the universe does 
not accelerate in long term future and has no event 
horizon, so this problem does not appear. Moreover, 
if one takes into account ordinary matter, then the 
universe at the asymptotically large time slowly de- 
celerates. But one can easily describe by this model a 
universe that presently accelerates and has iln — 0.7. 
The present value of w in this case is about —0.7, but 
the average value w ( [l^ ) in this case is —0.79, which 
looks pretty safe from the point of view of the present 
observational data. 

It has been observed in |2l] that for the positive ex- 
ponential potentials related to higher dimensional su- 
pergravities the typical value is A^ > 6. However this 
constraint may be decreased due to the presence of the 
FI terms, as shown in p2| . An interesting question to 
be explored is the M/string theory origin of potentials 
with A^ < 2 which could provide us with cosmological 
models of the presently accelerating universe. 

In our case the model with V = Ae^'^ appears as 
a special limiting case of the N — 2 model considered 
in the previous section. But a purely exponential V = 
Ae^'^ appears in some other theories as well. 

Indeed, a model of = 2, d = 4 gauged supergrav- 
ity with 3 vector multiplets, without hypermultiplets, 
described in |39] , has exactly such a potential: 



-1/2^ 



1 



-R 



1 



2 



(33) 



where the other few scalar fields, which are coordi- 



nates of 



5^(1,1) 



50(2,2) 



coset space are con- 



(7(1) SO(2)xSO(2) 

stant. This Lagrangian can be deduced from the com- 
plete action given in . The relevant symplectic sec- 
tions are X^,Fa = va^SX^; X^X^t^^^ = 0,tjas = 
(1, 1,-1,-1). The gravity-dilaton part of the action 
is 



-R 



(2ImS)2 2ImS 



(34) 



The complete gauged supergravity action is consistent 
under condition that Im S < which provides a pos- 
itivity condition for the kinetic terms of scalars and 
vectors. For ImS = — c^^"^ the action (34) is reduced 
to ( p3| ) with a positive definite potential. 

The potential is defined in eqs. (9.54) and (9.58) 
of [p9|. The source of such positive potential is due 



to Fayet-Iliopoulos term with (P4 )^ = This is 
possible for an Abelian gauge group. For the positive 
sign the Abelian gauge group with FI terms must be 
neither the graviphoton, nor the gravidilaton, but one 
of the remaining vectors in the direction where 7744 = 
-1. 

The action (33h can be also identified with the 
model of Rcf. under condition that there is no 
gauging of the 50(2, 1) group, cq = 0, the rotation 
angle is vanishing, 6 = 0, and only FI terms with ei 
are present. In such case this model is reduced to the 
model proposed in and described by Eqs. (p3|), 
(25). The contribution to the potential surviving the 
limit eo — > and — *■ is coming from the term 
gV20(^gx/20c.^2^) _^ ^e^4>_ Since ei is a FI term 
1^ this clarifies the origin of the positive exponential 
potential in A^ = 2 supergravities. 



VI. A = 1 SUPERGRAVITY 

A. Supergravity quintessence model 

Now we take lessons learned in extended N > 2 
supergravities and use them as a guide for A^ = 1 
models which may give us a successful description of 
the recent past of our universe, including the acceler- 
ation period, and lead to particular prediction about 
the future. 

The first example is the one already known as a 'su- 
pergravity quintessence' model p4| with the potential: 



V 



e 2 — ■ 



(35) 



This model has a dS minimum at (/)cr = \foL- Inter- 
estingly, the mass of the scalar field near dS minimum 
is exactly the same as in the N = 2 model discussed 
above: 



m, 



(36) 



We solved the Friedmann equations for the cases 
with a = 4,6,11. In case of initial values of the 
fields not far from the dS attractor value, all features 
of the solution (a(t), r2£)(<), i(;(i)) are practically the 
same as in our N = 2 future dS model. However, if 
one considers the model with a = 11 |Q and take 
the initial value of the field far away from the mini- 
mum of the effective potential at = %/Tl, which is 
necessary to approach the tracker solution, then w{t) 
rapidly increases from —0.8 and for a long time stays 
above —0.2. The average value H) (|lj) in this case is 
about —0.65. This makes this scenario rather vulner- 
able El. 
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FIG. 7. The values of w for 0o = 1 in the supergravity 
quintessence model of Ref. as a function of redshift 
z. The value of w at present is about —0.8, in agreement 
with . But at jz > 1 the value of w is very high. 

A complete supergravity model, underlying the po- 
tential (|3^) is rather complicated and fine-tuned. In 
addition to the quintessence field, it involves several 
different fields charged under an abelian group, X- 
fields of charge 1, Y fields of charge -2. Various ad- 
ditional assumptions about the non-minimal Kahler 
potential and superpotential are necessary to finally 
come up with the everywhere positive supergravity 
quintessence potential (^. 

One could argue that this fine-tuning of a super- 
gravity model is a reasonable price for the tracker na- 
ture of the solution. However, when one modifies the 



inverse power law potential 



by multiplying it 



by an e~ as in Eq. (35), the potential acquires a 
minimum. All of the solutions at the time when Q jj 
becomes equal to 0.7 enter the immediate vicinity of 
this minimum. The value of the potential in this min- 
imum needs to be fine-tuned if we want this moment 
to coincide with the present time. 

Thus, this model requires one fine tuning, exactly as 
all other models with de Sitter minimum or with the 
exponential potential considered in our paper. In this 
sense, tracking in the model of does not seem to 
provide additional benefits because it does not resolve 
the fine-tuning problem. That is why we would like to 
look for simpler models based on iV = 1 supergravity 
that may also describe dark energy. 



B. Generic N = 1 supergravity and Polony i model 

The general potential of = 1 supergravity con- 
sists of an _F-term and a D-term. Here for simplicity 
we will restrict ourselves to the study of the simplest 
models without a D-term. In this case 

V{z) = e'^ [-3WW* + {V'W)g-\^ (VjW*)] . (37) 

Here JC is the Kahler potential and W{zi) is a su- 
perpotential of a chiral superfield z^. The covariant 



derivative on W is V'W = d'W + (d'JC) W . As a 
simplest example, consider the model with the mini- 
mal Kahler potential /C — zz* and the linear superpo- 
tential W{z) = iJ,'^{z + f3). This is the famous Polonyi 
model, which serves as a standard part of models with 
gravity mediated SUSY breaking, see e.g. ||2^. It is 
convenient to represent the complex field z as a sum 
of two canonically normalized fields, z = {(f> + ix)/V^- 
One can show that the minimum of the effective po- 
tential occurs at X = 0, so we can restrict ourselves to 
investigation of V^(0), where 



V2Mp 



(38) 



Here we temporarily restore Mp in our equations. 
Note that in the limit Mp —^ oo corresponding to 
global SUSY this potential becomes exactly ffat, i.e. 
the effective mass of the moduli field (j) vanishes. For 
P = 0{1) this potential has a minimum at = 0{Mp), 
and the mass m of the field (f> in the minimum is 

oifiyMp). 

Usually it is assumed that the Polonyi field leads 
to supersymmetry breaking in the observable sector. 
It is supposed to have mass m ~ 10'^ GeV. Since this 
mass appears only due to SUSY breaking because of 
the gravitational effects, it is protected with respect 
to radiative corrections. As one can easily see from 
Eq. (H), the value of H'^ ~ V{cj))/M^ for < 1 

and 101 <C V^Mp) is of the same order of magnitude 
as m'^{(j)) = V" . This is a rather general (though not 
unavoidable) property of moduli fields in A^ = 1 super- 
gravity: Typically one has ^ V{(j))/M^ - Iw^ (</>)!, 
see [E2| and references therein. This is very similar 
to the situation in extended supergravity described 
above. However, in N = 1 supergravity this relation 
is not rigidly fixed. In particular, one can fine-tune the 
parameter /3 to be equal to 2 — In this case the 
effective potential vanishes in its minimum, which is 
necessary to avoid having an enormously large cosmo- 
logical constant A ~ and still have supersymmetry 
breaking with 1113/2 ^ m ^ 10^ GeV. 

However, if one can achieve supersymmetry break- 
ing with vanishing cosmological constant by a different 
method, then instead of the field with mass m ^ 10'^ 
GeV (or in addition to it) one can consider an ultra- 
light Polonyi field with m ~ 10"^"^ eV as a part 
of the dark energy hidden sector. Then the abso- 
lute value of the effective potential in its minimum 
for P = 0(1) is 0(10"^^") ~ po, so there is no reason 
to fine-tune /3 to be 2 - \/3 « 0.268, see Fig. |. It is 
important, that for V{(t)) > and /3 = 0(1) one has a 
generic relation ~ m?{4>), just like in extended su- 
pergravity. As a result, many of our results obtained 
in application to the N = 8 and N = 2 models can be 
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easily extended to a large class of = 1 supergravity 
models. 



V 




FIG. 8. Polonyi field potential for /3 = (symmetric 
potential with a minimum at V > 0), /3 = 0.2 (minimum 
at F > 0), /3 = 2 - ^/3 (minimum a.t V = 0), and P = 0.4 
(minimum at V < 0). In all of these cases \V"\ = 0(iV|) 
for 101 ^ 1- The potential is shown in units of m^; the 
field is shown in units of Mp. Potentials with /3 < can 
be obtained by the change (j> — > —(p. 

In particular, Polonyi model with < 2 — \/3 w 
0.268 leads to asymptotically de Sitter universe, just 
like in = 2 model of Ref. Ig^l • The fine-tuned model 
with = 2 — Vs asymptotically leads to Minkowski 
space. Meanwhile all models with 2 — \/3 < \(3\ lead 
to a collapsing universe. However, just like the = 8 
model, the = 1 models with 2 - VS < \(3\ < 0.5 
can describe dark energy in an accelerating universe, 
see Figs. ^ - These figures show the results of 
calculations where we for definiteness took the initial 
value of the field (p equal to — —1. 




t 



FIG. 10. The future of the universe in the Polonyi 
model for (5 — OA and (j>o ~ — 1- The present time cor- 
responds to t — 0. The universe in this regime would ac- 
celerate for the next 20 billion years, but then eventually 
collapse. In the models with /3 > 0.4 the collapse occurs 
much earlier. 




FIG. 11. The values otQo for the Polonyi field potential 
as a function of redshift z. The lower (horizontal) curve 
corresponds to /3 = 0, the next one - to /3 = 0.2, then 
(3 = 2- V3, and /3 = 0.4. 




FIG. 9. The scale factor a(t) for the Polonyi field po- 
tential. At large t the upper curve corresponds to l3 — 0, 
the next one - to (3 = 0.2, then = 2 - \/3, and = 0.4. 



w 




FIG. 12. The values of w for the Polonyi field potential 
as a function of redshift z. The lower (horizontal) curve 
corresponds to P = 0, the next one - to /3 = 0.2, then 
/3 = 2 - ^/3, and /3 = 0.4. 



Thus, even the simplest models based on A^ = 1 
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supergravity with m ~ 10"'^^ eV can provide a natu- 
ral description of dark energy because of the generic 
relation |m^| ^ that is satisfied in many of such 
models. 



C. The axion dark energy model 

One of the popular models of quintessence based 
on N=l supergravity is the axion model. The original 
axion model of quintessence was proposed in |26|] . It 
has the potential: 

y = A[cos(0//) + C] . (39) 

In and in most of the consequent studies of the ac- 
celerated universe in this model it was assumed that 
C = 1. This means that the cosmological constant 
was assumed to vanish in the minimum of the po- 
tential, just as in many other quintessence models. 
The positive definiteness of the potential with C = 1 
and the fact that it has a minimum aX V = Q could 
be motivated, in particular, by global supersymmetry 
arguments. In supergravity and M/string theory the 
argument from global supersymmetry requiring C — 1 
is not longer valid. The derivation of the value of the 
constant term C in the axion potential from any fun- 
damental theory is not available. 

In p7[ the axion model of quintessence was stud- 
ied using the arguments from M/string theory. These 
arguments were based on the properties of the mem- 
brane instantons. Some assumptions were made con- 
cerning the supcrpotcntial and Kahler potential in 
the effective = 1 supergravity potential ( |37| ) mo- 
tivated by M/string theory. The axion potential was 
presented there in the form 

y^Acos ((/.//) , (40) 

without any constant part. Perhaps it is possible to 
add to this potential a positive constant C, but here 
we will analyse this model without making any such 
additions. This potential has a maximum at = 0, 
V{{)) = A, and a minimum at = /tt, F(/7r) = —A. 
Therefore the universe collapses when the field (f) rolls 
down from the top of this potential. The axion model 
with the potential (|o|) given in |Q has the properties 
that are qualitatively similar to those described in our 
paper if one assumes that / = 0{Mp). Indeed, the 
curvature of the effective potential in its maximum at 
= is given by 

^"--f.= -j^^o ■ (41) 
If one takes, for definiteness, / = Mp = 1, one finds 
= -A = -3iJo , (42) 



and for / = Mp/V2 one has 

TO^ = -2A = -6H^ , (43) 

exactly as in the A^ = 8 model. 

It is quite possible to describe dark energy using this 
model. The results of the investigation of this model 
are very similar to our results for N — 8 supergravity, 
so we will not reproduce here the figures for a{t), 
and w. Instead we present in Figure |3| the range of all 
possible values of A and (jjo which allow the universe 
(described by the model with / = Mp) to live longer 
than 14 billion years (the area under the thick (red) 
curve). A is given in units of po and (j) in units of 
Mp. Note that this graph is periodic and symmetric 
under refiection 0o ~^ ~4'o- Therefore it is sufficient 
to present our results for < (po < /"'• 
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FIG. 13. The region below the thick (red) line contains 
all possible A and (j>o corresponding to the total lifetime of 
the universe greater than 14 billion years as predicted by 
the axion model (|^. The (blue) dashed line A ^ 3.6po 
separates this region into two equal area parts. The region 
below the thin (black) curve corresponds to all universes 
with the lifetime greater than 28 billion years, i.e. to the 
universes that would live longer than 14 billion years after 
the present moment. The area under this curve is 4 times 
smaller than the area under the thin (red) curve. 

To conclude, the study of extended supergravities 
attracted our attention to the properties of other re- 
alistic models, like the axion one. As we see from the 
Figure |l^, the set of parameters that allow the uni- 
verse to have the total lifetime of more than 28 billion 
years (i.e. to live another 14 billion years from now 
on) is 4 times smaller than the set of parameters that 
allow our universe to live only 14 billion years. This 
implies that the axion model with V = Acos(0//) 
and / = 0{Mp) can describe dark energy, but as in 
many previous cases, the universe can be expected to 
collapse within the next 10 - 20 billion years. The 
future collapse allows the anthropic considerations to 
be applied to such models, which in turn leads to a 
reasonable order of magnitudeprediction of the value 
of the cosmological constant . 
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VII. DISCUSSION 

In this paper we studied various possibiUties to de- 
scribe dark energy in supergravities and the future 
of the universe in such models. Most of the previous 
works on this subject concentrated on the attempts to 
obtain inverse power law tracker potentials in iV = 1 
supergravity p^j2^ . These models required introduc- 
tion of rather complicated superpotentials and Kahler 
potentials of several superfields and many assump- 
tions. 

In this paper we have chosen a different strategy. 
First of all, we considered several toy models based 
on extended supergravity which may have closer re- 
lation the M-theory. In addition we considered the 
simplest dark energy models in iV = 1 supergravity. 
All of these models share a very interesting feature: 
The absolute value of the effective mass squared of the 
scalar field responsible for the dark energy of the uni- 
verse is of the same order as the effective potential of 
this field, which implies that |m^| = 0{H^). Whereas 
in the phenomenological models of quintessence this 
property usually is required for their consistency, in 
supergravity this feature is rather common and some- 
times it is even unavoidable. In particular, in all 
known versions of extended supergravity with de Sit- 
ter solutions with the Hubble constant H, scalar field 
masses are always quantized: — nH^, where n are 
some integers which can be either positive or negative 
|l^j2^,|l^. For example, in all versions of iV = 8 su- 
pergravity of this type one always has a tachyon field 
with mass m? = = — 677g = — 2Vo, where Vb is the 
value of the effective potential in its extremum corre- 
sponding to de Sitter solution, and is the curvature 
of the potential at that point |l^Jl^. This makes all 
models with \m?\ — 0{V) interesting candidates for 
the role of the dark energy. 

In all models with de Sitter vacuum state with ~ 
Vq > the field slowly rolls to the minimum of the 
effective potential and the universe eventually reaches 
de Sitter state with = Vb/3. 

On the other hand, all iV 8, iV = 4, iV = 2 



models with 



— V" — 

— Vq — 



'QHn = — 2Vn, as well as 



= 1 Polonyi models with (3 > 2 — -s/S, describ- 
ing the present state with « 0.7, lead to the 
following generic prediction concerning the future of 
the universe: Our universe is going to collapse within 
the time comparable to its present age t lA billion 
years. Similar results are valid for the axion model of 
quintessence with the potential V = A[cos (0//) -I- C] 
for C < 1. 

The possibility that a flat universe may collapse 
in the future was known for quite a while |^ - ||l7|] . 
However, this possibility seemed to be rather extrava- 
gant, especially in view of the observations suggesting 
that the universe is accelerating. And even though 
we knew that the universe may collapse in a distant 



future, there was no reason to expect that this may 
happen relatively soon. 

Now the situation becomes quite different. Among 
all models of dark energy based on extended super- 
gravity only the N = 2 model of Ref. |20| leads to 
a stable de Sitter space in the future; all other mod- 
els lead to a collapse. Among the simplest A^ = 1 
Polonyi models only the models with f3 < 0.268 lead 
to a stable de Sitter space. Similarly, among all axion 
models with the potential V = A[cos {<p/ f) + C] only 
the models with C > 1 lead to a stable dS space. All 
other models predict that our universe should collapse 
within the next 0(10) billion years. 

Of course, all of the models considered in our pa- 
per are just toy models. We assumed that the dark 
energy hidden sector can be successfully incorporated 
into the theory of elementary particles and that the 
cosmological constant problem in the observable sec- 
tor can somehow be solved. But this is a general issue 
with all models of dark energy. On the other hand, the 
new class of models may provide an unusual solution 
to the coincidence problem. 

Indeed, the total lifetime of the universe in A^ = 
8 theories with de Sitter solutions |jl^,|l3|,|l8| is 
0{H-^) - 0(A-i/2). This lifetime can be few times 
greater, but only if the initial value (/)o of the field is 
exponentially close to = 0. Thus, large values of A 
are forbidden since they do not allow the long-living 
universes. 

If one considers all combinations of (j)Q and A com- 
patible with the total lifetime of the universe ttot > 14 
billion years and assume that all such combinations 
are equally probable, one finds that the value of A 
should very close to po ^ 10"^^", and the most prob- 
able time before the future collapse of the universe is 
0(10) billion years. This may provide a solution to 
the fine-tuning/coincidcncc problem for A and f2u si- 
multaneously predicting the typical time-scale of the 
global collapse of the universe in the models based on 
N — 8 supergravity | |3^ . Similar conclusion is valid 
for the simplest N — 1 Polonyi models with (3 > 0.268. 

Interestingly, all models that lead to the "dooms- 
day prediction" have some features that may allow us 
either to rule them out or to estimate the time remain- 
ing until the global collapse. For example, the N = 8 
model with (po < 0.3 and f2_D — 0.7 would lead to a col- 
lapse within the next 10 billion years. Meanwhile, the 
model with (j)o = 0.35 would lead to a collapse within 
the next 7 billion years, see Fig. |l|. However, this 
model leads to a maximal value of flo ~ 0.65, which 
is possible but not particularly favoured by the recent 
observations The models with (j)o > 0.4 would 
place the doomsday even much closer to the present 
moment, but they lead to Qjj < 0.56, which is incom- 
patible with the present cosmological measurements 
oi ilu. 

An additional information can be obtained by the 
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measurements of w. In all models predicting the 
doomsday in the near future, the value of w tends to 
grow significantly at small z, i.e. at the present epoch, 
see Figs. ||, ^ It is difficult to determine the time 
dependence of w using CMB experiments alone ||^, 
but one can combine them with the supernova obser- 
vations, counts of galaxies and of clusters of galaxies, 
and with investigation of weak gravitational lensing 
|4^,^,Q. Even in this case it is difficult to find the 
equation of state of dark energy w{z). However, if one 
takes the cosmological models based on supergravity 
seriously and realizes that our future is at stake, one 
gets an additional strong incentive to develop obser- 
vational and theoretical cosmology. It was never easy 
to look into the future, but it is possible to do so and 
we should not miss our chance. 
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